A test of the spin-orbit sum rule for actinides using an ab-initio calculation 
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The expectation value of the angular part of the spin-orbit-coupling (SOC) operator {I ■ s) and 
the branching ratio for the A'4,5 X-ray-absorption spectroscopy in 5-Pu, a-U and face-centered-cubic 
(fee) Th are calculated within the framework of the density-functional theory in combination with 
the local-spin-density approximation (LSD A), the generalized-gradient approximation (GGA) and 
the LDA+(/ method. The SOC contribution is calculated in terms of the first- and the second- 
variational schemes. A strong variation in the magnitudes of the calculated (Z • s) operator along the 
actinide series is found. The results show that the SOC sum rule for 5-Pu and a-U is reasonably 
valid, whereas there is a sign mismatch for the case of fee Th. The calculated branching ratios for 
the case without SOC in the valence shell strongly deviate from the statistical value. 

PACS numbers: 78.70.Dm; 71.27.-|-a; 71.15.Mb 
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The field of actinides attracts considerable inter- 
est in terms of experiments and theoretical investiga- 
tions. Recently, many unique phenomena of the ac- 
tinide series were discovered, for example, the Pu-based 
superconductivityii'^ or the phonon dispersion in 5- 
Pu^i^, which were ascribed to the nature of the 5/ 
electronic states. In order to explain these interesting 
properties, it is necessary to understand the electronic 
structure. It was demonstrated experimentally that the 
spin-orbit coupling (SOC) could not be neglected in the 
Hamiltonian for the 5/ states of the actinides^. 

Van der Laan and Thole^ proposed a probe for the 
SOC interaction, which is based on the branching ra- 
tio of the core-valence transition in X-ray absorption 
spectroscopy (XAS). The corresponding sum rule for the 
4d — > 5/ transition, which is relevant for actinides, reads 
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Here, {I ■ s) denotes the angular part of the SOC operator 
for a particular valence shell with Nh holes. The branch- 
ing ratio B for an electron transition from a d core shell 
is defined as: 
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where 15/2 and 73/2 are the total absorptions for the d^/2 
and d^/2 levels, respectively. The quantity Bq is the cor- 
responding branching ratio if there was no SOC in the 
valence shell. In the case of a negligible core-valence in- 
teraction the statistical value of Bq for a d core shell is 
3/5. There are several problems related to the appli- 
cation of the sum rule (^ on experimental data. The 
number of holes Nh for a particular shell is not exactly 
known, while the quantity Bq is not measurable, there- 
fore it is necessary to rely on some estimates or hints from 
calculations. Furthermore, the validity of Eq. (^ is not 
guaranteed since it is based on a free-atom model. For 
example, it was demonstrated^ that the SOC sum rule 



held well for the Ad and 5d metals but that it was violated 
for the 3d and 4/ metals due to the core- valence interac- 
tions. The aim of the present paper is to test the validity 
of the SOC sum rule for actinides by comparing the SOC 
expectation value {I ■ s) calculated directly from the elec- 
tronic structure with the value obtained from the calcu- 
lated branching ratio via Eq. (^ . A theoretical approach 
has several advantages over an experiment, because such 
a test cannot be performed on measured data, and both 
the experimentally uncertain quantities, Nh and -Bq can 
be easily determined with a calculation. 

The calculations for (5-Pu, a-U and the face-centered- 
cubic (fee) phase of Th were performed within the 
framework of the density functional theory by apply- 
ing the Wien97 codsiS, which adopts the full-potential 
linearized-augmented-plane-wave (FLAPW) methodii. 
The experimental values of the lattice parameters were 
used. The convergency tests implied that 1000 ((5-Pu and 
fee Th) and 864 (a-U) k points in the full Brillouin zone 
(BZ) were enough for a calculation of the corresponding 
quantities with the desired accuracy when the modified 
tetrahedron mcthod^S was used for the BZ integration. 
The plane- wave cut-off parameters were 6.7 Ry ((5-Pu), 
9.8 Ry (a-U) and 5.6 Ry (fee Th). 

Since the correlation effects are important for 
actinidesiii^ii^, the influence of the local-spin-density 
approximation (LSDA)i&, the generalized-gradient ap- 
proximation (GGA)ii and the LDA+C/ methodiSiiS was 
investigated. The parameters U and J, which appear in 
the LDA-|-[/ scheme, were set ioU — 2eV and J = 0.5eV 
for all three systems. This universal choice was based on 
values found in the literature for plutonium (Ref. Ilfj) 
and uranium (Ref. l20l) . and it was simply generalized to 
the case of thorium. The density-functional theory can 
yield a magnetic solution for actinides^^, hence the up 
and down spin states were distinguished. However, for 
all three systems the resulting magnetic moments were 
zero within the prescribed accuracy. The criterion for the 
self-consistency was the difference in the charge densities 
after the last two iterations being less than 10~'^e/(a.u.)'^. 

Special attention was paid to the calculation of the 



SOC contribution. The standard approach within the 
FLAPW method is to apply the second variational 
schem o^^i^? . In this method the eigenvalues and the 
eigenvectors are calculated in two steps within each it- 
eration. First, a scalar relativistic Hamiltonian (with- 
out the SOC term) is diagonalized. Then, the lowest 
of the resulting scalar-relativistic eigenvalues and eigen- 
vectors are used as a restricted basis set for diagonaliz- 
ing the full Hamiltonian with the included SOC term. 
While the second-variational method is sufhciently accu- 
rate, for example, for the Sd transition metals, it might 
break down for he avy elements like the actinides, as dis- 
cussed in Refs. 123124 A possible way to improve the re- 
sults is to extend the basis of the second-variational step 
by including relativistic pi/2 local orbitalg-^. An alterna- 
tive, more straight-forward, albeit less time-effective way, 
is to perform the calculations in terms of the so-called 
first-variational scheme, where the full Hamiltonian, in- 
cluding the SOC term, is diagonalized in a single step, 
using the full basis set of the linearized-augmented plane 
waves. In the present paper, the first- and the second- 
variational schemes were applied. The SOC term was 
simply set to zero for the calculation of the branching 
ratio Bq. 

The XAS spectra /i5/2(e) and /i3/2(e) as a function 
of the photon energy e were calculated using Fermi's 
golden rule in a nonrelativistic dipole approximation that 
is based on the evaluation of the matrix element for the 
operator p • e with e denoting the polarization vector of 
the light. The total absorptions n^/2 and nj/2, which are 
required for the calculation of the branching ratios B and 
Bq (2), were obtained from the integrals: 
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where E-p and Eq denote the Fermi energy and the upper 
edge of the 5/ valence band, respectively. The latter 
quantity is defined with the number of holes N^ and the 
corresponding density of states: n^f{e) as: 
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The results are presented in Table 1. Common to all 
the considered systems is a strong deviation of Bq from 
the statistical value of 3/5 ascribed to a non-negligible 
coupling between the core and valence electrons. The 
calculated Bq is almost a constant for S-Pu, a-\J and fee 
Th, regardless of the method applied for the calculation 
of the exchange-correlation potential, although there is a 
slight increasing trend from plutonium to thorium. This 
might be due to the mixing of the 4^5/2 and 4^3/2 core 
states, because the splitting between the two core lev- 
els grows along the actinide series, as demonstrated in 
Fig. 1. The magnitudes of the expectation values of 
the SOC operator (/ • s) obtained from the calculated 
branching ratios B and Bq via Eq. (1) in general over- 
estimate the corresponding quantities calculated directly 
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FIG. 1: The calculated total XAS spectra with (solid lines) 
and without (dashed lines) SOC in the 5/ valence shell for 
5-Pu (upper part), a-U (middle part) and fee Th (lower part) 
obtained by applying the LSDA and the second-variational 
method. 



from the electronic structure. The agreement is the best 
in the case of S-Pu, within ~ 4% to ~ 27%, depending 
on the calculational details. While there is almost no dif- 
ference between the results obtained using the LSDA or 
GGA, and the second-variational scheme yields almost 
identical results as the diagonalization of the full Hamil- 
tonian with the included SOC term, the LDA+U values 
differ considerably. The magnitudes of the {I ■ s) opera- 
tor obtained by applying the first-variational method are 
larger than the LSDA or GGA values. This is the only 
case where the magnitude of the directly calculated {I ■ s) 
operator is smaller than the branching-ratio value, but 
the agreement is very good, within 4%. However, the 
situation is the opposite for the case where the second- 
variational scheme was used for the calculation of the 
SOC term. The magnitudes are smaller than those ob- 
tained using the LSDA or GGA, and the magnitude of 
the branching-ratio SOC expectation value exceeds the 
directly calculated quantity by about 27%. 

The expectation values of the SOC operator for a-U 
are one order of magnitude larger than the correspond- 
ing quantities for (5-Pu, as it would be expected on the 
basis of the XAS-spectra plot in Fig. 1, where the differ- 
ence between the B and -Bo curves is substantially less 
pronounced when compared to the case of (5-Pu. The 
magnitudes of the branching-ratio values are about two 
times as large as those of the directly calculated quan- 
tities. The difference in the results obtained using the 
LSDA, GGA or the LDA+U method is subtle, as is the 
influence of choice of the method for including the SOC 
term. 

The expectation value of {I ■ s) for fee Th is reduced 
by another order of magnitude in comparison with S-Pu 
and a-U. While the directly calculated quantity remains 
negative, the values obtained by applying Eq. (I) are 
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S-Pu 


LSDA 


-0.548 


(-0.555) 


-0.606 


(-0.600) 


8.785 


(8.802) 


0.698 




GGA 


-0.547 


(-0.554) 


-0.605 


(-0.599) 


8.778 


(8.796) 


0.698 




LDA+U 


-0.672 


(-0.310) 


-0.645 


(-0.393) 


8.729 


(8.853) 


0.699 


a-\J 


LSDA 


-0.035 


(-0.035) 


-0.066 


(-0.059) 


11.430 


(11.514) 


0.701 




GGA 


-0.034 


(-0.035) 


-0.068 


(-0.066) 


11.518 


(11.496) 


0.700 




LDA+17 


-0.033 


(-0.038) 


-0.066 


(-0.066) 


11.537 


(11.619) 


0.699 


fccTh 


LSDA 


-0.002 


(-0.002) 


0.006 


(0.008) 


13.512 


(13.524) 


0.702 




GGA 


-0.002 


(-0.002) 


0.005 


(0.008) 


13.516 


(13.528) 


0.702 




LDA+U 


-0.001 


(-0.001) 


0.006 


(0.009) 


13.612 


(13.622) 


0.702 



TABLE I: A comparison between the directly calculated expectation value of the SOC operator (w^^") per number of holes N^ 
in the 5/ shell and the quantity — | (S — Sq), obtained from the calculated XAS spectra. The quantities B and So represent 
the branching ratios for the cases with and without SOC in the 5/ valence shell, respectively. The values in brackets were 
obtained by applying the second-variational scheme, while the rest of the values resulted from the first-variational treatment 
of the SOC term. 



positive, so that the SOC sum rule fails in this case. The 
reason is the absorption without SOC is larger than the 
absorption with SOC at the A^4 edge for fee Th, in con- 
trast to the case of 6-Pu or a-U, as can be seen in Fig. 
1. If the statistical value of 3/5 for Bq was used instead 
of the result of the calculation, the sign of {I ■ s) from 
the sum rule would be negative, as it is the directly cal- 
culated quantity, but the order of magnitude would be 
wrong. The details of the calculation, namely the choice 
of the exchange-correlation potential and the SOC treat- 
ment, do not have a significant influence on the results 
for fee Th. 

In conclusion, the applicability of the SOC sum rule for 
actinides was investigated theoretically by comparing the 
quantities obtained from the calculated branching ratios 
with the quantities calculated directly from the electronic 
structure. The agreement is the best for 6-Pu, while the 
application of the sum rule yields the wrong sign for fee 
Th. It has to be noted that the corresponding quantities 



for (5-Pu are two orders of magnitude larger than those 
for fee Th, hence the failure of the sum rule for the latter 
case might be partly ascribed also to the limited accuracy 
of the calculation. It was found that neither the choice of 
the exchange-correlation potential (the LSDA or GGA) 
nor the presence of the Coulomb repulsion (the LDA+U 
method) had a substantial influence on the quantities 
under consideration. The second-variational scheme was 
proven to be adequate for the calculation of the SOC con- 
tribution for most of the considered cases. The exception 
is the LDA+U calculation on 6-Pu. The problem here is 
with the mixed, spin-up and spin-down, matrix elements 
of the orbital-dependent potential, which are supposedly 
too large to be taken into account by means of the second- 
variational scheme. 

The experimentalists should bear in mind that the 
branching ratio Bq for the case where the SOC interac- 
tion in the valence band is switched off strongly deviates 
from the statistical value. 
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